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SWIM-LIKE MOTION OF BODIES IMMERSED IN AN IDEAL FLUID 

Marta Zoppello 1 and Franco Cardin 2 


Abstract. The connection between swimming and control theory is attracting increasing attention 
in the recent literature. Starting from an idea of Alberto Bressan [8] we study the system of a planar 
body whose position and shape are described by a finite number of parameters, and is immersed in a 
2-dimensional ideal and incompressible fluid in terms of gauge field on the space of shapes. We focus 
on a class of deformations measure preserving which are diffeomeorphisms whose existence is ensured 
by the Riemann Mapping Theorem. We face a crucial problem: the presence of possible non vanishing 
initial impulse. If the body starts with zero initial impulse we recover the results present in literature 
(Marsden, Munnier and oths). If instead the body starts with an initial impulse different from zero, 
the swimmer can self-propel in almost any direction if it can undergo shape changes without any bound 
on their velocity. This interesting observation, together with the analysis of the controllability of this 
system, seems innovative. 
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Introduction 


In this work we are interested in studying the self-propulsion of a deformable body in a fluid. This kind of 
systems is attracting an increasing interest in recent literature. Many authors focus on two different type of 
fluids. Some of them consider swimming at micro scale in a Stokes fluid HBiiaiH , because in this regime 
the inertial terms can be neglected and the hydrodynamic equations are linear. Others are interested in bodies 
immersed in an ideal incompressible fluid HPFPF1 and also in this case the hydrodynamic equations turn 
out to be linear. 

We deal with the last case, in particular we study a deformable body -typically a swimmer or a fish- immersed in 
an ideal and irrotational fluid. This special case has an interesting geometric nature and there is an attractive 
mathematical framework for it. We exploit this intrinsically geometrical structure of the problem inspired 


by 39 40 and 32 , in which they interpret the system in terms of gauge field on the space of shapes. The 


choice of taking into account the inertia can apparently lead to a more complex system, but neglecting the 
viscosity the hydrodynamic equations are still linear, and this fact makes the system more manageable. The 


same fluid regime and existence of solutions of these hydrodynamic equations has been studied in 18 regarding 
the motion of rigid bodies. 

We start from an early idea of Alberto Bressan [8] and some unpublished developments, according to which the 
shape changes can be described by a finite number of parameters. These kind of systems, where the controls are 
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precisely given by further degrees of freedom of the systems, have been first studied deeply by Aldo Bressan, 
see e.g. [9}|IT]. In this framework we show that the composed system “fluid-swimmer” is Lagrangian geodesic. 
Next, coupling this fact with some techniques developed in 33 , we are able to show that the kinetic energy 

and introduced in control 


38 


of the system (i.e. the Lagrangian) is bundle-like, a concept by Bruce Reinhart 

theory by Franco Rampazzo in [3T . This leads us to express the equations of motion as linear control equations, 
where any quadratic term is vanishing, radically simplifying our final analysis on the system. The geometric 


construction of the control dynamic equations follows substantially the line of thought of 13 29 


At a first glance, the deformations of the swimmer are naturally given by diffeomorphisms, that are infinite 
dimensional objects. By considering a planar setting and making use of complex analysis, as suggested in Mp3 


the Riemann Mapping Theorem plays a crucial role in describing the shape changes of the swimmer. It turns 
out that the diffeomorphisms can be parametrized by appropriate complex converging series. In the literature 
other authors exploit the same way of describing the shape changes by conformal maps, for example in 
in the environment of the Stokes approximate regime or in 


33 


14 


in the case of an ideal and irrotational fluid, 
in which they take into account only a finite number of terms to represent the diffeomorphisms. We follow 
substantially an analogous approach to merge this idea with the setting of Alberto Bressan. The choice of using 
a finite number of parameters means that the kind of deformations that we consider is more restricted but still 
enough to describe a wide range of swimmers. 

In order to have a more manageable system that the one in 


33 , we establish a connection between the use of 


complex and real shape parameters. We show that, if we consider small shape changes, a well precise choice 
of the real and imaginary part of the shape parameters leads to obtain exactly the same deformation proposed 
which use a rather different parametrization governed by suitable small deformation. Therefore we get 


32 


a description of our system with a finite number of parameters/coordinates, which is useful to apply the idea of 
controlling the shape coordinates to steer the swimmer between two different configurations. In this environment 
we recover the well known Scallop Theorem j3] in the case in which we suppose to have only one real shape 
parameter. Thanks to the idea of using a finite number of parameters we can reduce our dynamic equations 
to a control system. The controllability issue has been recently linked to the problem of swimming [2][3|[22|[2 7 [ 
since it helps in solving effectively motion planning or optimal control problems. 

We point out that in the Stokes regime there are interesting results, for example in 


27 


the controllability of a swimmer which performs small deformations around the sphere, or in 23 


the authors study 
in which he 


considers the swimming mechanism as a ‘broken-line’-like structure, formed by an ordered sequence of finitely 
many sets. Even in [2 16 they study the controllability of a slender swimmer composed by N links immersed 
in a viscous fluid at low Reynolds number. 

In the present work we deal with the ideal and irrotational case, neglecting viscosity. In particular differently 
from what was done in other works, we focus our attention on a crucial problem: the presence of an initial 
impulse. The case of zero initial linear momentum is studied in literature concerning systems of different nature: 
both in the multi-particle or many-bodies field, 19-21 , and shape changing bodies, 30 32 -34 , as the equation 


of motion are a driftless affine control system whose controllability can be studied using classical techniques. 
Instead, the case of a non vanishing initial impulse leads us to a more complex system since the equations of 
motion involve also a non zero drift term and their controllability is more tricky to study. Therefore we have 
two contributions to the motion of the system: the first one that is purely geometrical and determined by the 
structure of the problem, and the second one, strictly linked to the presence of a non vanishing initial impulse. 
The controllability of this kind of systems is studied in detail, and among other facts it is worth noting that we 
need at least three real shape parameters to make the system controllable. 

We have three state parameters, three conjugate variables and at least three controls. Despite the evident 
complexity of the computations linked to this number of variables, we managed to obtain interesting results. 


The plan of the paper is the following. In Section 1 we present in some detail the geometric aspects useful 
to formulate our problem. The proper geometrical setting of the swimmer in a 2-dimensional fluid is faced in 
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Section 2. Section 3 contains an exhaustive study, in a complex setting, of the deformation of the body, together 
with the construction of the equation of motion. We deal with all the controllability issues in Section 4. 


1. Preliminaries 

This section covers some auxiliary mathematical topics, in particular from Lie groups, fiber bundles and 
connections that we shall need later. This summary will be helpful to set the notation, fill in some gaps, and 
to provide a guide to the literature for needed background. 

1.1. Lie Groups 

Let us start from some geometric and algebraic notions on Lie groups, that arise in discussing conservation 
laws for mechanical and control systems and in the analysis of systems with some underlying symmetry. 

Definition 1.1. A Lie group is a smooth manifold G that is a group with identity element e = gg^ 1 = g~ l g, 
and for which the group operations of multiplication, ( g,h ) i-a gh for g,h £ G, and inversion, g i-a g~ x , are 
smooth. 

Before giving a brief description of some of the theory of Lie groups we mention an important example: the 
group of linear isomorphisms of M™ to itself. This is a Lie group of dimension n 2 called the general linear group 
and denoted by GL(n,R). The conditions for a Lie group are easily checked. This is a manifold, since it is an 
open subset of the linear space of all linear maps of K™ to itself; the group operations are smooth, since they 
are algebraic operations on the matrix entries. 

Definition 1.2. A matrix Lie group is a set of invertible nxn matrices that is closed under matrix multiplication 
and that is a submanifold ofM nxn . 

Lie groups are frequently studied in conjunction with Lie algebras, which are associated with the tangent 
spaces of Lie groups as we now describe. 

Definition 1.3. For any pair of n x n matrices A, B we define the matrix Lie bracket [A, B ] = AB — BA. 

Proposition 1.1. The matrix Lie bracket operation has the following two properties: 

(i) For any nxn matrices A and B, [B, A] = — [A, B] (skew-symmetry). 

(ii) For any nxn matrices A, B, and C, 

[[A, B\,C\ + [[ B , C], A] + [[C, A], B] = 0. (This is known as the Jacobi identity.) 

As is known, properties (i) and (ii) above are often thought as the definition of more general Lie brackets 
(than AB — BA) on vector spaces called Lie algebras. 

Definition 1.4. A (matrix) Lie algebra g is a set of n x n matrices that is a vector space with respect to the 
usual operations of matrix addition and multiplication by real numbers (scalars) and that is closed under the 
matrix Lie bracket operation [■,■]. 

Proposition 1.2. For any matrix Lie group G, the tangent space at the identity TjG is a Lie algebra. 

As usual, for matrix Lie groups one denotes e = I 

We now describe an example that plays an important role in mechanics and control. 

The plane Euclidean Group 
Consider the Lie group of all 3 x 3 matrices of the form 


R d 
0 1 


( 1 . 1 ) 
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where R £ SO( 2) and d £ R 2 . This group is usually denoted by SE(2) and is called the special Euclidean 
group. The corresponding Lie algebra, se(2), is three-dimensional and is spanned by 


Ai — 



/0 0 V 
A 2 = 0 0 0 

\0 0 0, 


'0 o o> 
A 3 = I 0 0 1 
,0 0 0 ; 


( 1 . 2 ) 


The special Euclidean group is of central interest in mechanics since it describes the set of rigid motions and 
coordinate transformations on the plane. Let G be a matrix Lie group and let g = TjG be the corresponding 
Lie algebra. The dimensions of the differentiable manifold G and the vector space g are of course the same, and 
there must be a one-to-one local correspondence between a neighborhood of 0 in g and a neighborhood of the 
identity element I in G. An explicit local correspondence is provided by the exponential mapping exp : g K > G, 
which we now describe. For any A £ R nx " (the space ofnxn matrices). exp(A) is defined by 


xp(A) := I + A+ ^A 2 + ^A 3 + ... 


(1.3) 


This map for SE(2) can be defined by the exponential of the elements of the Lie algebra se(2). More precisely 


exp(dAi) = 


exp(a’A2) = 


Since [Ai, Aj\ = 0 for all i, j = 1, 2,3, we have that V(0, x, y) € R 3 = g = se(2) 



(1.4) 


(1.5) 


exp(Mi + xA 2 + yA 3 ) = exp(6b4i) exp(au42) exp(yAs) 


that is clearly elements of SE( 2). 

We now define the action of a Lie group G on a manifold Q. Roughly speaking, a group action is a group 
of transformations of Q indexed by elements of the group G and whose composition in Q is compatible with 
group multiplication in G. 

Definition 1.5. Let Q be a manifold and let G be a Lie group. A left action of a Lie group G on Q is a smooth 
mapping $ : G i —> Q such that 

(i) 4>(e, q) = q for all q£Q, 

(ii) 4>((jr, $(/i, q )) = <&((//)., q) for all g,h £ G and q £ Q, 

(iii) 4>(g, •) is a diffeomorphism for each g £ G. 

A Lie group acts on its tangent bundle by the tangent map. We can consider the left or the right action of 
G on g by: T e L g f or T e R g £, where L g and R g denote left and right translations, respectively; so if g = g(t) is 
a curve in G, then there exists a time dependent £(•) eg such that 


9(t) = T e L g{t) f{t) = g(t)£(t) 


( 1 . 6 ) 


and similarly for the right action. 

Given left action of a Lie group G on Q, 4 1 : G x Q —x Q, and f an element of the Lie algebra g then 
<lA : R x Q —> Q : ( t,q) i—> 4>(exp tf, q) is a flow on Q , the corresponding vector field on Q is called infinitesimal 
generator of <f> corresponding to £, is denoted by f,Q(q) 

£ q ( q ) = ^$(expf^,g)| t=0 . 


( 1 . 7 ) 
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1.2. Fiber Bundles and Connections 

Fiber bundles provide a basic geometric structure for the understanding of many mechanical and control 
problems. 

A fiber bundle essentially consists of a given space (the base) together with another space (the fiber) attached 
at each point, plus some compatibility conditions. More formally, we have the following: 

Definition 1.6. Let S be a differentiable base manifold and G a Lie group. A differentiable manifold Q is 
called principal fiber bundle if the following conditions are satisfied: 

1 G acts on Q to the left, freely and differentiably: 

$:GxQ->Q (1.8) 

writing <F(g, q) = • q = g ■ q. 

2 S = Q/G and the canonical projection n : Q —> S is differentiable 

3 Q is locally trivial, namely every point s £ S has a neighborhood U such that 7r -1 (t/) C Q is isomorphic 
to U x G, in the sense that q £ 7r _1 ({/) 1 —> (n(q), (f>(q)) £ U x G is a diffeomorphism such that <f> : 
7r _1 (t/) —X G satisfies <f>(g ■ q) = g<f>(q),\/g £ G 

An important additional structure on a bundle is a connection. Suppose we have a bundle and consider 
(locally) a section of this bundle, i.e., a choice of a point in the fiber over each point in the base. We call such 
a choice a “field”. The idea is to single out fields that are “constant”. For vector fields on a linear space, for 
example, it is clear what we want such fields to be; for vector fields on a manifold or an arbitrary bundle, we have 
to specify this notion. Such fields are called “horizontal”. A connection is used to single out horizontal fields, 
more precisely fields which live in a subspace of the the tangent space, and is chosen to have other desirable 
properties, such as linearity. 

Definition 1.7. Let (Q,S,n,G) be a principal fiber bundle, the kernel ofT q n denoted by V q := {v £ T q Q\T q ir(v) = 0}, 
is the subspace ofT q Q tangent to the fiber through q and is called vertical subspace. A connection on the prin¬ 
cipal fiber bundle is a choice of a tangent subspace H q C T q Q at each point q £ Q called horizontal subspace 
such that: 

(1) T q Q = H q ®V q 

(2) For all g £ G and q £ Q, T q <h g ■ H q = H g . q 

(3) H q depends differentiably on q 


Hence, for any q £ Q, we have that T q w determines an isomorphism H q = T^( q) S: for all T q Q 9 v = vy q +Vh q 
and we have that T n ( q )(v) = Vn q £ S. In other words the choice of an horizontal subspace can be seen also as 
the choice of a vector valued “connection one form” which vanishes on the horizontal vectors. 

It follows the definition 


Definition 1.8. An Ehresmann connection A is a vector valued one form such that 

(i) A is vertical valued: A q : T q —> V q is a linear map for each point q £ Q 

(ii) A is a projection: A(v) = v for all v £ V q . 


In the special case in which ( Q,S,tt,G ) is a principal fiber bundle the previous conditions on A : TQ —> g 
read: 

(i) A(£q(< 7 )) = £ for all £ £ g and q £ Q, where £g(g) is the infinitesimal generator of the left action of G 

on Q ([1J|. ' 

(ii) A is equivariant: 


A(T q ^ g (v))) = Ad g (A(v)) 
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for all V £ T q Q and g £ G where <f> g denotes the given action of G on Q and where Ad denotes the 
adjoint action of G on 0 defined as 

Ad g .- T e (^Lg O Rg~ l) . 0 t 0 . 

Therefore it is evident that the horizontal subspace H q is the kernel of A q . 
n the case in which there is a metric h{q) in our manifold Q , we have a special way to define the horizontal 
subspace: it is the orthogonal with respect to the metric to the vertical subspace. 

Hq = {w £ TqQ : (w,h(q)v) = 0,Vu £ V q } (1.9) 

In this special case our connection A is called mechanical connection (see j3Tj and therein references). We now 
would like to express the connection in coordinates, in order to do this we first introduce the following definition 

Definition 1.9. Let us consider the following diagram 


Q 



S D u 


The functions like a are sections and we call T(U, Q) the set of all sections from U in Q. 

Alternatively often a connection is introduced as a derivation V as follows. Let V be a map 

V : T(Q) T(Q ® T*S) such that 
V(fTi + 02 ) = V(cTl) + V((7 2 ) 

V( fa ) = /V(ct) + a 0 df if / is a C°° function. 

Let now e be a local basis of sections of the principal fiber bundle, in this basis the connection one-form A can 
be expressed as 

e a A% = Veb a, b = 1 ■ • • dim(Q). 

If we change basis in r(Q), say e = ef l, the connection A changes, i.e. 

eA =Ve = V(efi _:1 ) = (Ve)fT 1 + edfT 1 = eAfT 1 + edfT 1 
= efMfW 1 + ef IdLl^ 1 

therefore A and A satisfy the following relation 

A = flAfl- 1 + fldLl- 1 ( 1 . 10 ) 

Let u(t ) be a smooth curve in S passing through the point P = m( 0). Let q £ Qp = 7r -1 (P) be any point in 
the fiber of Q over P. We would like to find a smooth curve 7 (f) in Q such that 71(7 (t)) = u(t), 7 ( 0 ) = q , and 
7 '(t) £ (i.e., the tangent vectors to the curve 7 (t) are horizontal). 

From the usual theory of differential equations it follows that such a curve 7 (t) exists and is unique, at least 
locally at any point q £ Q (i.e., for small values of t). The curve 7 is called a horizontal lift of u. Regarding 
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the tangent vectors, for any q £ Q and any vector u £ T^r qj S there exists a unique vector v £ H q C T q Q such 
that T q TT : v iu. The vector v is called the horizontal lift of u. 

Given an Ehresmann connection we can define the horizontal lift of curves in <S, hence we can also define a 
notion of parallel transport that allows us to identify different fibers of Q. 



Figure 1. Horizontal lift of the curve u(t) 


Note that, in general, the parallel transport will be path-dependent. If we start with two different curves 
u\(t) and U 2 (t), such that tti(O) = 112 ( 0 ) = P and u\[f) = U 2 (i) = S, then the horizontal lifts 71 and 72 at a 
point q £ Qp will not meet, in general, at a point in the fiber Qs, he., we will have 71 (f) 7 ^ 72 (t). This gap 
on the fiber is called holonomy and depends on the choice of the connection and on the topology of the base 


manifold. In particular if it is connected the holonomy depends on the basepoint only up to conjugation 24 


2. Geometrical setting 

In this section we present the geometrical framework underlying dynamical control systems. We derive the 
equations of motion and discuss how to use the geometrical tools introduced before to gain informations on our 
system. 

2.1. Geometry of control equations 

In this subsection we derive the local dynamic equations for the control system ( Q , h , T) where J- is a smooth 
fc-dimensional foliation on Q , and h is the Riemaniann metric on the manifold Q, as done in 29 . As is well 
known, on a set U C Q adapted for the foliation, T coincides with the model foliation of K" by fc-dimensional 
hyperplanes. Let (j) '■ U —> R ra , <f>(P) = ( x , y ) be a local chart of Q in U, distinguished for T, so that cj) maps 
P\u into the trivial fibration 7 r(x,y) = y. Set q = (a :,y) £ Q; given a path u(t) £ 7 r(</>([/)), we suppose that for 
every t , the reaction forces that implement the (ideal) constraint y = u{t ) are workless with respect to the set 
V q{t) U = ker T q ( t )ir of the virtual displacements compatible with the constraint y = u(t). 

Let (Q, h, JF) be a foliated Riemanian manifold, let U C Q be an open set adapted for T an let q = (x, y). If 
T(q,q) = |() t /i(( 7 )q is the kinetic energy of the unconstrained system (Q.hjP), then the kinetic energy of the 
system subject to the time dependent constraint y = u(t) is T(x,u(t),x,ii(t)). The related dynamic equations 
are, in Lagrangian formalism 

d_ dT _dT_ lt 
dt dx dx 


( 2 . 1 ) 
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These can be put in Hamiltonian form by performing a partial legendre transformation on the x- variables. 
When we identify y with u(t) and y with J \u{t ), the above Lagrange equations are equivalent to 


x = 


dH 

dp 


( X,p,U,il ) 


dH 

P = U ’ U ) ■ 


( 2 . 2 ) 


We call these equations control equations. Let 


= xfCx + x t My + + y l By 


(2.3) 


be the local block representation of the metric h in 4>{U), where C,B are symmetric and invertible respectively 
k x k and (n — k) x (n — k ) matrices. 

To every q £ U denote with H q U the subspace orthogonal to V q U = kerT 9 7r with respect to h. Referring to 
the local expression of h in U, it is easy to see that H q U is the space orthogonal to the vectors (e*, 0), = i...„ with 
respect to the metric h . 


H q ll = {(±, y) € T q U such that C(q)x + M(q)y = 0} . 


Therefore H q U can be equivalently assigned through the V q [/-valued connection one form defined in 1.8 


d 


A(q) = (dx + C(q)dy) (g) -A- where (see 2.6 ) C = C 


(2.4) 


whose kernel and range are respectively H q U and V q U. Now we consider the orthogonal splitting of a vector 
into its horizontal ad vertical components 

v = v v + v h — A(q)v + hor(T q nv) = (x + Cy , 0) + (— Cy , y) 

Using the above decomposition, we get the induced splitting of the kinetic energy metric tensor into its vertical 
and horizontal part: 

h(q)dq ® dq = C(q)A(q) 0 A(q ) + K{q)dy 0 dy (2.5) 

where K{q) = B — 

Definition 2.1. The Riemannian metric h is bundle-like for the foliation J- iff on a neighborhood U with 
adapted coordinates (x 7 y) the above orthogonal splitting of g holds with K = K{y). 

The importance of this notion will be clear in the following subsection ( 2.1.1) . Using this notation we want 
to rewrite the control equations. 

From 

dT „ 

P = vjY = Cx + My 

we obtain 

x = C~ 1 p — C~ 1 Mu = C~ 1 p — Cu 
Substituting (2.6 1 in (2.31 and recalling that — ^ = §y we have 


( 2 . 6 ) 


dH dT 1 t dC~ 1 t dC 

P = -frc = dx=-2 P -^ P + P ^x U ' 


1 ,d(B-M t C~ 1 M) . 
2 U dx u 


(2.7) 


Therefore the control equations are 


| x = C 1 p — C l Mu 

[P — -?P flrr P + P U ^ “ 


( 2 . 8 ) 
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We now introduce, following 29 
of a surjective submersion 7r 
with pq : T*Q 


the global version of the above dynamic equations when Q is the total space 
Q —> S. Let VQ be the vertical subbundle and V*Q the dual of VQ. Denote 
Q the cotangent projection and set n := tt o pq, tt : V*Q —> S. If (x,y) are local fibered 


coordinates on Q, (x,y,p) are local fibered coordinates on V*Q. Moreover, denote with z = (x,p) the local 
coordinates on the ff-fiber over y. Now, to every y £ S, n~ l (y) is a fiber canonically simplettomophic to 
T* (7r 1 (?/)), representing the phase space of the constrained system restricted to the 7r-fiber over y. 


T*Q 


V*Q TQ 

( x,y,p ) (x,y,x,y) 



y 


Suppose that a control vector field u is given on S and that the path u(t) is an integral curve of u. Then the 
dynamic equations (2.61 and (2.7) are the local expression of a vector field Du over V*Q that projects on u by 
tt. Moreover the field Du is tangent to the fiber of tt only if the control is trivial: u vanishing. Let us suppose 
that the control is given by a curve u : [ti, < 2 ] —> S in S that is the integral curve of the vector field u. Thus the 
movement of the system is described by a differentiable curve 7 : [ti, ^ 2 ] —> Q such that 71 ( 7 ( 4 )) = u(t). Note 
that (S : [ti,t 2 ] —> TQ is the natural increase of the curve 7 in the fiber tangent to Q. Composing 4* with 
the Legendre transform CtQ —> T*Q and with the projection r : T*Q —> V*Q we obtain the parametric curve 
7 = r o C o ^ : [ti, t 2 \ —> V*Q which represent the evolution of the system taking into account the control. 

Let horQ : TS — > TQ denote the horizontal lift of the Ehresmann connection, introduced in the previous 
section, and pq the cotangent projection, using the above definitions we introduce the function 


Ku ■ V*Q —¥ R Ku o pq(q) = ( horQ(q){u)) t h{q)hor Q {q){u) 

Theorem 2.1. To every control vector field u on S, the corresponding dynamic vector field Du can be expressed 
as the sum of three terms: 


Du = X Ho - X Kii + hor(u) 

(2.9) 

v d 1 t dc - 1 d 

Xh - =C p 0x 2 P dDdv 

(2.10) 

„ 1 .,sk . a 

~ Xk - = 2 “ & 

(2.11) 

, ... / d d t dC d . . 

/,or(«) = (--C- +P —-) u 

(2.12) 


with 
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where Xh 0 is the Hamiltonian vector field, corresponding to the case of locked control, X]is the Hamiltonian 
vector field on V*Q associated to Kf and hor is the horizontal lift of an Ehresmann connection on 7 r : V*Q —> S 
entirely determined by 7 r and the metric. These equations are exactly the control equations (2.8 1 . 


Proof: 13 


□ 


2.1.1. The importance of initial impulse 

In what follows let us suppose that the metric h is bundle like. 


Proposition 2.2. The control system (2.8 1 is of two different types depending on the value of the initial value 
of the x conjugate variables p. 


(1) Case p(0) = 0 


The system (2.8) is an affine non linear driftless control system; 
(2) Case p(0) 7 ^ 0 


The system (2.8) is an affine non linear control system with drift. 
Proof: Since we have supposed to have a bundle like metric we have that 

d K{y) d{B - M t C~ 1 M)(y) 


dx 


dx 


= 0 . 


Therefore the control equation (2.8 1 becomes 


| x = C l p — C 1 Mu 

\p = -hp t ^?p+p t ^ 


(2.13) 


Case p(0) = 0. 

The function p(t) = 0 is the unique solution of (2.13 >2 according to the Cauchy theorem. Thus (2.13) 1 becomes 
a driftless control system. 

x = —C~ 1 Mu 


It is clear that this last equation is entirely determined by the connection (see (2.4)). Therefore in the case of 
null initial impulse case only the geometry of the system determines its motion. 

Case p(0) / 0. 


In this case the equation (2.13 >2 has no trivial solution that is p(t) 7 ^ 0. Thus (2.13) 1 is a non linear control 


system with drift determined exactly by the presence of a non zero p 

x = C~ x p — C^ 1 Mii 


The presence of the drift is crucial because in this case the motion of the system is determined both by the 
connection (given by the geometry) and by the impulse, that is non zero. This proves the importance of the 
initial value of p. □ 


In this work we analyze both the cases. The one with zero initial impulse is well studied in literature for 

The one with p(0) 7 ^ 0 is becoming of increasing interest since the presence of the 


many systems 19-21 1132 33 


impulse influences deeply the motion, as we have seen. We deal with this problem that is more complex and 
tricky to study because of the presence of the drift. 
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2.2. Geometric and dynamic phase 

As we have seen, in the general theory, connections are associated with bundle mappings, which project 
larger spaces onto smaller ones. The larger space is the bundle, and the smaller space is the base. Directions 
in the larger space that project to a point are vertical directions. The connection is a specification of a set 
of directions, the horizontal directions, at each point, which complements the space of vertical directions. In 
general, we can expect that for a horizontal motion in the bundle corresponding to a cyclic motion in the base, 
the vertical motion will undergo a shift, called a phase shift, between the beginning and the end of its path. The 
magnitude of the shift will depend on the curvature of the connection and the area that is enclosed by the path 
in the base space: it is exactly the holonomy. This shift in the vertical element is often given by an element of 
a group, such as a rotation or translation group, and is called also the geometric phase. Referring to what 
said in the previous subsection, the motion is determined only by the geometrical properties of the system if 
it starts with zero initial impulse. In many examples, the base space is the control space in the sense that the 
path in the base space can be chosen by suitable control inputs to the system, i.e. changes in internal shape. In 
the locomotion setting, the base space describes the internal shape of the object, and cyclic paths in the shape 
space correspond to the movements that lead to translational and rotational motion of the body. 

Nevertheless the shape changes are not the only ones to determine a net motion of the body. More generally, 
this motion can always be decomposed into two components: the geometric phase, determined by the shape 
of the path and the area enclosed by it, and the dynamic phase, driven by the internal kinetic energy of 
the system characterized by the impulse. It is important to stress the difference between the two phases. The 
geometric phase is due entirely to the geometric structure of the system. Instead the dynamic phase is present 
if and only if the system has non zero initial impulse or if the impulse is not a conserved quantity, in our context 
we refer to what is explained in subsection |2.1.1| More precisely if the curvature of the connection is null, not 
necessarily the system does not move after a cyclic motion in the base: a net motion can result if the system 
starts with non zero initial impulse, and this motion is entirely due to the dynamic phase. 



Figure 2. Geometric phase and dynamic phase. 

Figure ([2]) shows a schematic representation of this decomposition for general rigid body motion. In this 
figure the sphere represents the base space, with a loop in the shape space shown as a circular path on the sphere. 
The closed circle above the sphere represents the fiber of this bundle attached to the indicated point. Given 
any path in the base (shape) space, there is an associated path, called the horizontal lift, that is independent 
of the time parametrization of the path and of the initial vertical position of the system. Following the lifted 
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path along a loop in the shape space leads to a net change in vertical position along the fiber. This net change 
is just the geometric phase. On top of that, but decoupled from it, there is the motion of the system driven 
by the impulse, (if it is not zero) which leads to the dynamic phase. Combining these two provides the actual 
trajectory of the system. 

2.3. Gauge potential 

Let us consider a planar body immersed in a 2 dimensional fluid, which moves changing its shape. For the 
moment we do not specify the kind of fluid in which it is immersed that can be either ideal and incompressible or 
a viscous one with low Reynolds number. Our aim is to show that the motion of this deformable body through 
the fluid is completely determined by the geometry of the sequence of shapes that the idealized swimmer 
assumes, and to determine it. This idea was introduced by Shapere and Wilczek in 40 39 and developed 


14 , where they apply geometrical tools to describe the motion of a deformable body in a fluid, focusing their 


attention on the Stokes regime. 

The configuration space of a deformable body is the space of all possible shapes. We should distinguish between 
the space of shapes located somewhere in the plane and the more abstract space of unlocated shapes. The latter 
space can be obtained from the space cum locations by making the quotient with the group of rigid motions 
in the plane, i.e declaring two shapes with different centers of mass and orientation to be equivalent. The first 
problem we wish to solve can be stated as follows: what is the net rotation and translation which results when a 
deformable body undergoes a given sequence of unoriented shapes? The problem is intuitively well posed: when 
a body changes its shape in some way a net rotation and translation is induced. For example, if the system is 
composed simply by the body, its net rigid motion can be computed by making use of the law of conservation 
of momentum, if instead the body is immersed in an ideal incompressible fluid this motion can be found by 
solving the Euler equations for the fluid flow with boundary condition on the surface of the body with the shape 
corresponding to the given deformation. 

These remarks may seem straightforward, but we encounter a crucial ambiguity trying to formulate the problem 
more specifically. Namely how can we specify the net motion of an object which is continuously changing shape? 
To quantify this motion it is necessary to attach a reference frame to each unlocated shape. This is equivalent 
to choosing a standard location for each shape; more precisely to each unlocated shape there now corresponds a 
unique located shape. Once a choice of standard locations for shapes has been made, then we shall say that the 
rigid motion required to move between two different configurations is the displacement and rotation necessary 
to align their centers and axes. In what follows we shall develop a formalism, already used in 40 39 , which 


ensures us that the choice of axes for the unlocated shapes is completely arbitrary and that the rigid motion on 
the physical space is independent from this choice. This will be clear soon below. 

For a given sequence of unlocated shapes So(t), the corresponding sequence of located shapes S(t) are related 
by 

S(t) = K(t)S 0 (t) (2.14) 

where TZ is a rigid motion. This relation expresses how to recover the located shapes S(t) given the unlocated 
ones, i.e. So(t). It is clear that we are dealing with a fiber bundle: the located shapes S(t) live on the big 
manifold Q = SE( 2) x S and the unlocated ones, So(t), live on the base manifold obtained by the quotient of 
the manifold Q by the plane euclidean group SE( 2), i.e S = Q/SE( 2). 

To make ( 2.14| ) more explicit we introduce a matrix representation for the group of Euclidean motions, of 
which TZ is a member. A two dimensional rigid motion consisting of a rotation R followed by a translation d 
may be represented as a 3 x 3 matrix 

r R d> 


[R,d} = 


0 1 


(2.15) 


where R is an ordinary 2x2 rotation matrix, d is a 2 component column vector. This is the matrix representation 
of the plane euclidean group action SE( 2) on the manifold Q where the located shapes S(t) live on. 

Now in considering the problem of self propulsion we shall assume that our swimmer has control over its form 
but cannot exert net forces and torques on itself. A swimming stroke is therefore specified by a time-dependent 
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sequence of forms, or equivalently unlocated shapes So(/). The located shape will then be expressed exactly by 
formula ( 2.14| ). 

Our problem of determining the net rigid motion of the swimmer thus resolves itself into the computation of 
7 Z(t) given So(t). In computing this displacement it is most convenient to begin with infinitesimal motions and 
to build up finite motions by integrating. So let us define the infinitesimal motion A(t) by 


dlZ , , dlZ . . 

= iziiz - 1 —) = iza 


dt 


dt 


(2.16) 


In this formula we can recognize the differential equation corresponding to formula (1.6), from which we under¬ 
stand that A take values in the Lie algebra of the plane euclidean group: g = se(2). For any given infinitesimal 
change of shape A , formula (2.161, describes the net overall translation and rotation which results. We can 
integrate it to obtain 

rt-2 


lZ(t 2 ) = lZ(t 1 )Pexp / A(t)dt 
L Jt ! 

where P denotes a reverse path ordering, known in literature as chronological series [lj: 


(2.17) 


/'exp 


A(t) dt 


= 1 + 




A(t) dt - 


A(t)A(t ) dt.dt 


ti<t <t<t 2 


The assignment of center and axes can be arbitrary, so we should expect that physical results are independent 
of this assignment. How does this show up in our formalism? A change in the choice of centers and axes can 
equally well be thought of as a change (rigid motion) of the standard shapes, let us write 


s 0 = n(s 0 )s 0 


(2.18) 


The located shapes S(t) being unchanged, (2.141 requires us to define 40 39 

K(t) = K(t)n~ 1 {S 0 {t)) 


From this, the transformation law of A follow 


A = flAfl -1 + n 




dt 


(2.19) 


from which we can recognize the transformation laws (1.10) of an Ehresmann connection called also Gauge 


potential. Our freedom in choosing the assignment of axes shows up as a freedom of gauge choice on the space 
of standard shapes. Accordingly the final relationship between physical shapes is manifestly independent of 
such choices. 

Our aim will be to compute this gauge potential A £ se( 2) in function of the unlocated shapes So that our 
swimmer is able to control. 


3. Swimming in an ideal fluid 

We focus on a swimmer immersed in an ideal and incompressible fluid. The dynamical problem of its self 
propulsion has been reduced to the calculation of the gauge potential A. In our model we assume that the 
allowed motions, involving the same sequence of forms will include additional time-dependent rigid displace¬ 
ments. In other words the actual motion will be the composition of the given motion sequence So(t) and rigid 
displacements. 
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3.1. System of coordinates 

Let (0,ei,e2) be a reference Galilean frame by which we identify the physical space to R 2 . At any time 
the swimmer occupies an open smooth connected domain B and we denote by T = R 2 \ B the open connected 
domain of the surrounding fluid. The coordinates in (0,ei,e2) are denoted with x = {x\,X 2 ) T and are usually 
called spatial coordinates. Let us call (— X 2 ,x±) T = ar 1 . 

Attached to the swimmer, we define also a moving frame (0*,eJ,e2). Its choice is made such that its origin 
coincides at any time with the center of mass of the body. This frame represents the choice of the axes in the 
space of unlocatecl shapes. As we have shown before, the computation of the net rigid motion of the swimmer 
due to shape changes is independent from this choice that accordingly is arbitrary. The fact that this frame has 
always its origin in the center of mass is a matter of convenience: indeed this choice, and others (see Remark 

We denote by x* = (x^xlf) 1 the related so called body coordinates. In this frame and at any time the swimmer 
occupies a region B* and the fluid the domain JF* := R 2 \ B*. 

We define also the computational space , that is the Argand-Gauss plane which we will need only to perform 
explicit calculations, endowed with the frame (O, Ei, E 2 ) and in which the coordinates are denoted z = {z\, zf) T . 
In this space D is the unit disk and O := R 2 \ D. 

Remark 3.1 (Minimal Kinetic Energy). In order to simplify the calculations, since, as said in the previous 
section, the choice of the body frame is arbitrary, we use the one of the minimal kinetic energy. It is the one 
in which the velocity of the center of mass is null. This condition is clearly satisfied if its origin coincides with 
the center of mass at any time. Moreover in this frame the angular momentum with respect to the body frame 
has to be null. 

Remark 3.2. The orientation of the body frame remains arbitrary and does not effect the fact that it is the 
frame of minimal kinetic energy. One of the most used conventions to define a possible orientation of such a 
system is to choose as axes the eigenvectors of the moment of inertia of our body. Obviously as we have said 
before this choice does not effect the located shape, since it is independent on the choice of the frame. 


3.1 1 , tell us that the body frame is the one in which the kinetic energy of the body is minimal 26 


3.2. Shape changes 


Banach spaces of sequences. Inspired by 33 , we denote any complex sequence by c := (ck)k> 1 where for 


any k > 1, Ck '■= Ofc + ibk £ C, a*,, bk £ R. Most of the complex sequences we will consider live in the Banach 
space 

S := {( Ck)k>i : ^ K\a k \ + \b k \) < + 00 } 


fc>i 


endowed with its natural norm ||c||,g := + IM)- This space is continuously embedded in 


T := {(cfc)fc>i 


sup 

zed-D 


kckZ k | < + 00 } 

k> 1 


whose norm is ||c|| 7 - := sup zedD \ Y^k> 1 k°kZ k |, where D is the unit disk of the computational space. 

Definition 3.1. x We call T> the intersection of the unit ball ofT with the space S. 

This space will play an important role in the description of the shape changes that will follow. 
Finally we introduce also the Hilbert space 


ht {{ck)k> 1 : E K\ a k\ 2 + M 2 ) < + 00 } 

k> 1 


whose norm is ||c||y := \/J2k> 1 ^(l°fc| 2 + IM 2 )- According to Parseval’s identity we have 
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^ fc|cfc| 2 < fc 2 |c fc | 2 = — / kc k e~ lke ) 2 d9 < sup | ^ kc k z k 

k> 1 fc> 1 ^ fc>l zeaD fc> 1 

Therefore we have the following space inclusions 


S cTcU 

We have introduced these spaces because they will be crucial in the description of the shape changes of the 
idealized swimmer. 


3.2.1. Description of the shape changes 

Following the line of thoughts of [33] and 14 the shape changes of the swimmer are described with respect 
to the moving frame (O*, e^, e^) by a C 1 diffeomorphism x(c), depending on a shape variable c £ D which maps 
the closed unit disk D of the computational space onto the domain B* in the body frame. The diffeomorphisms 
x(c) allows us to associate to each sequence c a shape of the swimmer in the body frame. We can write, 
according to our notation, that for any cgP (see definition 3.11, 


x(c):CD5^1 2 E(0*,e(,e;) (3.1) 

and B* = x(c)(D). 

We now explain how to build the map x( c ) for an y given sequence c, see Fig [3] 



Figure 3. The physical space and the body frame. 

Theorem 3.1 (Riemann Mapping Theorem). Let K. be a simply connected open bounded subset of C with 
0 £ 1C. Then there exists an holomorphic isomorphism f : D —¥ K, with /(0) = 0. Any other isomorphisms 
with /(0) = 0 are of the form z <—> f(rz) with r £ dD a rotation. All functions f can be extended to an 
homeomorphism of D onto K. if and only if OK. is a Jordan curve. 
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Defining Coo = {CUoo}, if O = C oo\D, from the isomorphism / we have also an isomorphism from D to the 
exterior F*; we apply to F* the inversion p(z ) := \ obtaining the open simply connected G, we find another 
Riemann- isomorphism g : D —> G with g(0) = 0. Then we consider h = p o g = 1 : D —> F*. The function g 

is injective around zero, therefore g (0) 0, it follows that h has a pole of the first order in zero and therefore 

has a Laurent expansion 


h{z) 


1 

z 


g"( 0) 
2 


OO 


+ C-kZ k 


(3.2) 


We now have the area theorem 


36 : if a function like h is injective on the punctured disk then we have 


5>|c fc | 2 <l (3.3) 

k —1 

If we want an isomorphism of O on F* we take 0(c)(z) = h(~) 

OO 

0(c) (z) = Z + Y1 % ( 3 - 4 ) 

k =1 2 


We now suppose that the boundary of B* is a Jordan curve, i.e. simple closed curve in the plane, therefore 
the function 0(c) can be extended to homeomorphism on the boundary. Now 0(c) : O —> F* can be extended 
continuously to all Coo setting in the interior of D 

X( c )(z) := z+ ^2c k z k , (z £ D) (3.5) 

k> 1 

Since z=(on dD we deduce that the following map is continuous in C for all c £ V: 


<L(c)(z) 


X(c)(0) if z&D 

0(c) (z) if z G O = Coo \ D 


(3.6) 


Proposition 3.2. For all c £ T>, x(c) : D —> B* and 0(c) : O —> T* are both well defined and invertible. 
Further, x(c)|d is a C 1 diffeomorphism, 0(c) |q is a conformal mapping and < f ) (c) is an homeomorphism form 
C onto C. 


Proof: 


33 


□ 


Remark 3.3. Despite the generality of the Riemann Mapping Theorem , the way in which we decided to represent 
our diffeomorphism, lead us to some restrictions. Indeed in order to be sure that also x(c) is well defined -from 
proposition 3.2- we need to impose the restrictive condition c £ D, see (3.31, meaning that the shape variables 
have to be finitely bounded for both the norms of S and T■ To summarize we can say that to use the shape 
variable c £ V allows us to describe all of the bounded non-empty connected shapes of the body that are not too 
far from the unit ball. 


3.3. Rigid motions 

The overall motion of our body in the fluid is, as said before, the composition of its shape changes with a 
rigid motion. The shape changes have been described in the previous subsection and, as we will see, the Gauge 
potential A described at the beginning depends only on the shape variable c that is 


A = A{ c) 


(3.7) 
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this will be evident in the next sections. 

The net rigid motion is described by an element of the planar euclidean group as explained in subsection |2.3[ 
More precisely it is given by a translation d, which is the position of the center of mass, and a rotation R of an 
angle 9 , that gives the orientation of the moving frame ( 0 *,eJ,e 2 ) with respect to the physical one. 

Let the shape changes be frozen for a while and consider a physical point x attached to the body. Then 
there exists a smooth function t (d(t),9(t)) such that the point’s coordinates in (0,ei,e2) are given by 
x = R(9)xq + d. Next compute the time derivative expression (d, 9). We deduce that the Eulerian velocity 
of the point is v^x) = 9{x — d)- 1 + d. It can be also expressed in the moving frame (0*,eJ,e2) and reads 
vj = (lix*) 1 - + d* where 

00=^00 < 3 - s > 

where 7Z is an element of the euclidean group SE(2) of pure rotation. 

Remark 3.4. Notice that d* is not the time derivative of some d* but only a symbol to expresse the velocity d 
in the body frame. 

Let us return to the general case where the shape changes are taken into account. We deduce that the 
Eulerian velocity at a point x of B is 

v{x) = 9(x - d)- 1 + d + TZ(e)x(c)[x(c)~ 1 {R(0) T (x- d))] 

where the last term represent the velocity of deformation and is computed taking into account that x = 
R t (9)x* = R T (9)x(c)(z). When we express this velocity in the moving frame we get 

v*{x*) = ( 9x + d*) + *(c)(x(c)-V)) , (3.9) 


which is more compact and will be useful in what follows. 


3.4. Dynamics for ideal fluid 

In this section we use some well known ideas developed for example in some works of A. Bressan [1], which 
further simplify the system of our idealized swimmer. 

We assume that the shape changes of our swimmer can be described by a finite number of shape parameters, 
= (ci, • • • , c m ), thus we can call q = (qi, • • • , q m + 3 ) = (d, 9, Ci, • • • , c m ). This choice is widely spread 
18 32 , and implies that we focus only on a class of deformations which consist of 


i.e. c 

in recent literature as in 

particular shape changes that are sufficient to describe a wide range of swimmer behaviors. Let us call \ the 
diffeomorphism which describes the superimposition of the shape changes with a rigid motion, more precisely 
y(q) (z) := [i?(0), d] oy(c) (z) Assuming that there are no external forces, we wish to derive a system of equations 
describing the net motion of the body due to the shape changes and of the surrounding fluid expressed in the 
moving frame. Let N = m + 3 and 

1 N 

T{ q,q) = 2 E (3.10) 

i,j=1 

describe the kinetic energy of the body. For simplicity, we assume that the surrounding fluid has unit density. 
Calling v = v{x,t) its velocity at the point x, the kinetic energy of the surrounding fluid is given by 


K = 



(3.11) 


If the only active force is due to the scalar pressure p , the motion of the fluid is governed by the Euler equation 
for non-viscous, incompressible fluids: 


Vt + v ■ Vv = — Vp 


(3.12) 
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supplemented by the incompressibility condition 


div v = 0. 


In addition, we need a boundary condition 

/ N 


V 


3x(q)(X). _ x (q)(z) 

/ v o_ Z ~\ 


k =1 


dqk 


lx (q)(z)l 


= o 


(3.13) 


x (aKg) 


—z = n { x )i i x = x(q)(~)j 2 € 9D) denotes the unit outer normal to the set x(c)(D) = B at the point 

x, and is computed making the complex derivative of the function x(q)(z) expressed in polar coordinates that 
is 

■ ^(x(q)(e i<T )) _ ia x (q)(e ig ) _ _ x(q)(z) 

|9 CT (x(q)(e i<T)) | lx(q)(e lff )| lx'(q)(z)| 

which states that the velocity of the fluid has to be tangent to the surface of the body. To find the evolution of 
the coordinate q, we observe that 


d dT _ dT 
dt dq k dq k 


+ F k 


k = l---N 


(3.15) 


where T is the kinetic energy of the body and Fj are the components of the external pressure forces acting on 
the boundary of B. To determine these forces, we observe that, in connection with a small displacement of the 
q L coordinate, the work done by the pressure forces is 


SW = ~Sq k 




ldD\ lx'(q)(*)l 9q k j 


) da 


(3.16) 


The equation of motion are 


d dT 
dt Qqk 


dT 

dq k 


( z ) \ p(x(q)(z))J(q)(z) da 




(3.17) 


idD \ lx (q)001’ dqk 

We now show that, in the case of irrotational flow, the coupled system can be reduced to a finite dimensional 
impulsive Lagrangian system. It is well known (see 18 28 ) that the velocity field of the fluid can be determined 
by setting v = 'Slip and solving the Neumann problem in the exterior domain 


A ip = 0 x € T 

n ■ SJip = n ■ v(x) U=x( q )(z) x <E dB 
\ip\ — > 0 |a;| —> oo 


(3.18) 


where the boundary condition reads 


^ x(q)00 v^<9x(q)^. 

n ■ v(x) = ■ 2 ^ ~^z—\ z )qk 


'lx(q)WI 


dq k 


(3.19) 


Let us now consider the function <p( q) : K 2 \ D —> R 2 defined by the composition of <j>(q)(z) with the rigid 
motion [R, d], that clearly on the boundary of D coincides with the function x(q). From the linearity of (3.18l 
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the solution will be linear in q. 


N 


i>{z, q,q) = ^^7k(z>q)9fc 


k =1 


The motion of the fluid can be obtained by solving the ordinary differential equation 


precisely 


^( q )(~) = £r(*»q.q)L=^ ( vow 


v L (x, q, q) = X] = E q)9fcL=^(q)(*) 


(3.20) 


(3.21) 


(3.22) 


This has to be true for all curve 1 9 1 1 —> q (t), thus we have 

a<^(q) _ 

'-'■'I lx=0(q)(z) 


dqk V ~ y dx L 

We now prove that the term of the equations of motion relative to the pressure forces is a kinetic term 

<^(q)(z) 94>{q) 


(3.23) 


Fk = - / 7T7 


-(^)p(^(q)( 2 :))J(q)( 2 :)da- = 


OD |0'(q)(z)| d 1k 

f wT(^( q )( z ))p(^( q )( z )) J ( q )( z ) dcr = 

IdD \(j)(q){z)\dx L 


S aB nL(x) i£ rdx 


applying the divergence theorem to (3.24) 


Lt dl = L (p ^i +V7t ' Vp) * = 


Vy/c ■ Vpdx = 




x€lT 

=0 

f d(f>{q) 

lze«. 2 \D 9q k 


■ (~v, t -v ■ Vv)j(q)(z) dz = 


94>(q) d - 

—uJ(q)(z) dz = 


Jzes. 2 \D 9qk dt 

f ,ddj>( q) L d 2 4>(q) ? 

= - / hr(—H- «i) - a fc gfc]J(q)(g) dz 

Jzen 2 \D dt oq k oqjdq k 

where J(q)(z) is the determinant of the jacobian matrix of the function 0(c) (z). Let us define 


(3.24) 




xGJ 7 


i N _ . . 

' dx ~ 


h3 


then 


F k = ~(- 


(i. [ d FLdx-l 

\dt J X £jr Qqk Jx<zj= dq k / 


d\v\ 2 \ d dTf dTf 

dx = — - 1 


dt dq k 9q k 


In conclusion the system body+fluid is geodesic of Lagrangian 
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In what follows for simplicity we will express all the quantities in the moving frame (O*, e*, e^), denoting the 

* 

total kinetic energy in this frame as T■ We will now compute explicitly the Lagrangian. Let us start with 
the kinetic energy of the swimmer. Since we have chosen the body frame as the one of minimal kinetic energy, 
according to Konig theorem, there is a decoupling between the kinetic energy of the body due to its rigid motion 


and that due to its shape changes, recalling (3.101: 


* body 111 

T '■= -m|d*| 2 + -I(c)6 2 + - 


x( c )(x( c ) (£*)) 


dm* 


(3.26) 


where 1(c) is the moment of inertia of the body thought as rigid with frozen shape, and the last term being the 
kinetic energy of deformation. It can be computed as follows: 


x( c )(x( c ) V*)) dm* = / \x(c)(z)\~ dm 0 = 7rp 0 ^ 

Jd 


k —1 


|Cfc| 

k + 1 


where we used the formula (3.51 to compute the integral. Note that accordingly to remark (3.2) the kinetic 


energy of the body in the frame (O*, e*^) does not depend on the orientation of the frame but only on its 
angular velocity. 


3.4.1. Kinetic energy of the fluid 

Since we are interested on the effect of the shape changes of the swimmer on the fluid, in this subsection we 
will compute all the quantities in the body frame. As we have seen in subsection |2.3| we can recover the rigid 
motion of the swimmer due to its deformation, exploiting the Gauge potential. 

The kinetic energy of the fluid reads 

T := lf \u*\ 2 dm}= | Vr\ 2 dm* f (3.27) 

There u* = ^7 if* and if* is the solution of the Neumann problem 

r Aif* = o 

| n(s*) ■ = n(i*) • «(®*)|i*= x ( q )(z) x* G dB* (3.28) 

[ \if*\ -)> 0 |x*| -> oo 


which is the same Neumann problem (3.181 expressed in the body frame. Indeed since the Laplacian operator 
is invariant under rototranslations, the function if*(x*) = if([R(9), d](x)) is harmonic. 

We will use complex analysis to compute the potential function if*. We define the function £(z) := 
tf*(<f(c)(z)), (z £ O), where if* is the potential function defined in ( |3.28 1 expressed in the moving frame 
and recalling (|3.4|) <f(c) (z) is the conformal map from O = C \ D to the external domain T* . According to 


classical properties of conformal mappings, the function £ is harmonic in O and the following equality holds: 


1 

2 



|Vi/>*| 2 dm*f 


±J o \V£\ a dm 0 f 


(3.29) 


The main advantage of this substitution is that £ is defined in the fixed domain O , whereas if* was defined in 
T* depending on c. 

In the moving frame is now easier to compute explicitly the boundary condition of the Neumann problem. The 
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outer normal to dB* is, recalling (3.141 


, , *X , . , *x 0 (c)( 2 ) 

n(x ) := m(x ) + m 2 (x ) = -2 , 

10 (c)( 2 )! 

where </> (c)(c) is the complex derivative of Recalling the following identity 

= d n i>f{x*) (x* = <t>(c)(z)) 

10 (c)( 2 )| 


(3.30) 


(3.31) 


and taking into account the expression (3.91 of v*, we deduce that the Neumann boundary condition (3.19l 
reads 


d n £(z) = • n = -dl$i(z(/) (c)(z)) - d^zcf (c)(z)) - 0 Q(</>(c)(z)z<fi (c)(z)) 


-®(x(c )20 (c)( 2 )). 


(3.32) 


This equality leads us to introduce the functions £J(c)(j = 1,2,3) and f d (c) as being harmonic in O and 
satisfying the following Neumann boundary conditions: 


dnCi ( c )( 2 ) = - 3 ?( 2 </>' (c)( 2 )), (3.33) 

dn&( c )( z ) = - 3 ( 20 '(c)( 2 )), (3.34) 

dnt 3 (c)( 2 ) = -3(&(c)(z)z0'(c)(z)), (3.35) 

d n ^ d (c)(z) =-Si(x(c)z(/>'(c)(z)), (zedD). (3.36) 


In this way we spilt the harmonicity and the Neumann boundary conditions of the function £ into the same 
properties for the functions £J(c) (j = 1, 2, 3) and f d (c). 

Next we have for all c £ V 


(c)( 2 )) = ^ dk$t(z k+1 </)' (c) ( 2 )) + 6 fc 3(2 fe+ V' (c)(z)) 

fc=i 

Proposition 3.3 (Potential decomposition). According to the Kirchhoff law the following holds in the Sobolev 
space H 1 (0): 


£(C) = <foi r (c) + d* 2 &( c) + 0&(c) + (£ d (c),c), (3.37) 

From the linearity of this expression with respect to d*, 9, c and since the gradient function preserves the 
linearity, we deduce that the kinetic energy of the fluid is a quadratic function of d*, 6 , c. 

3.5. The Gauge potential and the equations of motion 

According to what proved in the preceding section the Lagrangian of our system is a quadratic form in 
(d*, 0 , c), therefore it can be written in blocks as follows 

T(d*,e,c) = i(V T , 0 )M r (c) + 2 (d* T , 0 )N(c)c + c T M d (c)c) (3.38) 


where M r , N and M d play the role of the matrices C A4 and B , introduced in the section 2.1.1[ respectively 
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Remark 3.5. It is worth noting that in the physical space the kinetic energy is 

T(d, 8 , c, d, 0, c) 

When it is expressed in the body frame instead it becomes 

T ( c , d*, 8 , c) 

This expression does not depend on d and 8 due to the symmetry of our model with respect to the position and 
orientation of the body in the fluid. 


As we have seen in the first section we are interested in determining the Gauge potential A associated to our 
system which is 


. ^ , dlZ 

A = K~ 1 — = 
dt 




0 0 



(3.39) 


Since all the matrices M r , N and depend only on c, the kinetic energy is independent from d and 8 and 
the metric that it defines is bundle like (see Definition 2.11. In the principal fiber bundle SE( 2) x S —» S, the 


Gauge potential A depends on the kinetic energy, through the equation of motion, therefore also A does not 
depend on the state variables. 

We now need to determine (d*,d? 2 .8). In order to do this we compute the Hamiltonian associated to the 
Lagrangian function performing a partial legendre transformation on the q* variables. 

Before passing to formal calculations we recall how to interpret the connection introduced before in the cotangent 
bundle setting following the steps presented in subsection 2.1.1 This construction was presented also in [l3l29 . 


Let VQ be the vertical subbundle and V*Q the dual of VQ. Denote with pg : T*Q —> Q the cotangent 
projection and set ir := tv opg, n : V*Q — > S. If (d, 0, c) are local fibered coordinates on Q, (d,0,c, p*) are 
local fibered coordinates on V*Q. Suppose that a control vector field c is given on S and that the path c (t) is 
an integral curve of c. Then the equation of motion are the local expression of a vector field over V*Q that 
projects on c by n. Moreover the field Dc is tangent to the fiber of 7r only if the control is vanishing. 

Recalling that q = ( 91 , • • • , q m + 3 ) = (d, 6, c) 


'dT\ f d* 

P = = M - (C) 


N(c)c 


which defines the translational and angular impulses of the system body plus fluid. From this we obtain 


L (c)p* — IVE ” 1 (c)M(c)c 


(3.40) 


This expression is very convenient to study the motion of the shape-changing body since it gives the velocity 
with respect to the shape variable. 

It is easy to recognize the terms of the sum in which the control equations are split according to Theorem 

I O : 


X H 0 = 


0 


Xk, = 0 


and 


hor{ c) = 


ir-i 


(c)N(c)c 
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Note thus that we are exactly in the case of a bundle-like metric Therefore the equation of motion regarding 


the state variables are are exactly the ones given by formula (3.401. 


To obtain the equation of motion regarding the conjugate variables, we follow the method explained in 
Lamb 25 and Munnier 33 : we introduce P and II, the translational and angular impulses, as well as L and 


A, the impulses relating to the deformations: 



= M r (c) 



P 


* 


P + L 

n +A 



We start from the Lagrange equations 


d dT dT 
dt dqi dqi 


i = 1,2,3 


(3.41) 


and recall that T (c,d*,0,c) = T(c, R(9)d*, 9, c). 
Therefore recalling that dgR(O) = R(^)R(d) 


d dT dT 

dt f)d ()d 

d dT dT 

~dt ~dQ ~ ~dd 


s(fF flW ) = s< p+L >-"< p+L > J - 

p +L y = 

4(n + A)-d*(P + L)- L 


from these equations we get 


Pi 
dt \p* 2 


-P2 

Pi 


= 0 


-Pl - dipt + dipt = 0 


Therefore the equation of motion in the body coordinates are 


(3.42) 


(3.43) 


d* 


L (c)p* — M r 1 (c)N(c)c 


< pt = 0p* 2 
P*2 = -dpi 

Pt = d* 2 pt - dipt 


(3.44) 


Notice that these equations are exactly the ones presented in 
and the conservation of the impulse. 


32 which describe the evolution of the state 


3.5.1. Simplification of the deformations 

Since we are interested in studying small deformation around a circular shape, see remark |3.3[ in order to 
describe it we are interested only in expressing the distance from the origin of each point on circle’s boundary 
in function of the shape parameters c = a + ib. Let us consider deformation described by m shape parameters. 







24 


TITLE WILL BE SET BY THE PUBLISHER 


According to formula (3.5) the deformation can be written as: 

m 

X(c)(X) = 2 + E CkZ k for z € D 


k =1 


therefore the modulus of a point on the boundary described in polar coordinates by z = e l<T is given by 

m m 

|x(c)(Y)| 2 ={z + J2 c kZ k )(z + ^c k z k ) = 
k=1 k =1 

m 

(l + ^2(ak cos((fc + 1 )a) + b k sin((fc + 1 )<t)) (3.45) 

fc=i 

m 

+ y; (fflfe + ib k )(ah - ibh )) 

h,k=l 

taking the square root and using the Taylor expansion around c = 0 which corresponds to the circular shape 
we obtain 

h 
2 


lx(c)(z)| = ( 1 + E(y cos (( fc + 1 ) cr ) + y sin (( /c + 1 ) cr ))) + E°( C *) ( 3 - 46 ) 

fc=l fe=l 

where we can neglect all the terms of order grater or equal than 2 supposing afc, bk small for all k, for example 
of order e. 

For example in the case m = 2 taking 

ai = 2esi b-[ = 0 02 = 2eS2 = 2es3 (3.47) 

we find exactly the formula for the swimmer deformation proposed by Mason and Burdick in ( [32] ). 

F(er, s) = [l + e(si cos(2ct) + S 2 cos(3cr) + S 3 sin(3cr))] (3.48) 

This will be the formula for the deformation that we will take from now on. 

Remark 3.6. The fact that we are neglecting the terms of order e 2 is not in contrast with what said in j 33^ . 
They consider that the points can rotate on the boundary at high frequency producing a macroscopic identical 
deformation, leading to different dynamics. Instead our type deformations (3.461 do not allow this kind of 
behavior since it is a radial one. 


it— 1 


(s)p* - M r 1 (s)N(s)s 


We can express the equation of motion (3.441 using the real parameters Sk as shape parameters instead of 
Cfc obtaining: 

'?) 

Pi = Op* 2 (3-49) 

p* = -dp\ 

pi = dfp'l - d\p* 2 

where the matrices M r and N have the same physical meaning of the matrices M r and N but are expressed 
using the real shape parameters s. 

From now on we focus only on shape transformations near the identity, like (3.481 so that we can use real shape 
parameters to describe the deformation of the system. 








TITLE WILL BE SET BY THE PUBLISHER 


25 


3.5.2. Curvature of the connection: geometric and dynamic phase 

In this subsection we deal with the problem of having a net motion performing cyclical shape changes. Look¬ 
ing at equations (3.491 1 is evident that there are two contributions: the one of M r (s)p* which involves the 
impulse and — M r (s)N(s)s which is entirely geometrical. 





First, let us suppose that the system starts with zero initial impulse, i.e. p*(0) = 0. With this as¬ 
sumption the last three equations of the system (3.491 have as unique solution the null one therefore, 
the first term of equation (3.49 q vanishes and the infinitesimal relationship between shape changes and 
body velocity is described by the local form of the connection computed above. Moreover we take into 
account the reconstruction relation (3.8), which links the state velocity expressed in the body frame 
with the one expressed in the physical frame, 


V 0 


t r (s)N(s)s = -~gAi(s)s l 


(3.50) 


where g is an element of the planar euclidean group SE{ 2). From these we recognize the expression of 
the Gauge potential (2.161. 

We would like to find a solution for this equations that will aid in designing or evaluating motions that 
arise from shape variations. Because SE{ 2) is a Lie group this solution will generally have the form 

g(t) = <?(0)e*« 

where z £ se( 2), the Lie algebra relative to SE( 2). An expansion for the Lie algebra valued function 
z(t) is given by the Campbell-Hausdorff formula 


1t=f- 


17 ^ 


1 


z — A + -[A, A] + - [[A, A], A] + — [A, [A, A]] + 


(3.51) 


A(t) = f A(t)s(t) dr 
Jo 


To obtain useful results in the spatial coordinates, examine the group displacement resulting from a 
periodic path a : [0, T] —> R m , such that a(0) = a(T). Taylor expand A* about a(0) and then regroup, 
simplify, apply integration by parts and use that the path is cyclic 


z(a) = ~F iJ (a( 0)) [ ds l dJ + \(F lhk - [A u F ok ]){a{ 0)) [ ds'd.A d, k ■ ■■ 


(3.52) 


where 

Ftj = Aj t i — Aij — [Aj, Aj] 

is called curvature of the connection. 

For proportionally small deformations, the displacement experienced during one deformation cycle is: 

gdisp = e z{a) k. exp(—^Fjj(a(0)) J ds l ds J ) (3.53) 

If the curvature F is not null this displacement gives us the so called geometric phase that is the 


statement of the well-known Ambrose-Singer theorem 12 
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p’/O 


Let now suppose that the system starts with an initial impulse which is non zero. Thus the last 
three equations of (3.491 are not trivial. First of all we need to integrate this equations, which in 
function of the deformation s and s take the form 


fp* = (M7 1 (s)p*) 3 P2 - (A{s)s) 3 P 2 

< P*2 = — (Mr 1 (s)p*) 3 p5; + (A(s)s) 3 pt (3.54) 

[p*3 = (M(T 1 (s)p*) 2 p* - (A(s)s) 2 p! - (M“ 1 (s)p*) i Pa + (A( 


these can be solved once the shape s is prescribed as a function of time, and as before we choose a 
periodic shape path a : [0,T] —> R m ,with a(0) = a(T). 

Let us now consider the equation of motion regarding the state variables. We have both the contribu¬ 
tions: the geometrical one, already studied in the case with zero impulse, and also the one depending 
on the impulse p*. 

9 = ff( M r ^P* - 4li(s)s l ) (3.55) 

As before the integration of this term along a{t) gives 


9dis P = g( 0 )e 2(t) 


where 


z = Z+-[Z,Z] + - [[Z, Z\, Z) + - [Z, IZ, Z]] +. 


(3.56) 


Z := 




(t)p*(t) — A(t)s(t) dr 


In order to see that gdisp is effectively the sum of two contribution let us focus on the third equation of 
( |T55l ). It is 

6 = (M~ 1 (s)p*) 3 — (A(s)s ) 3 


from this we can easly recognize two terms. The first one integrated along a is 

/ (M7 1 (a(T))p*(r)) 3 dr, 

Jo 


(3.57) 


which value depends strictly on the evolution of the impulse p* given by equations (3.49 ) 2 — 4 - The second 
term is the geometric contribution analyzed in the previous section which depends on the curvature of 
the connection. 

Once we have integrated this system and obtained the time evolution of 9 we can solve also the ODEs 
regarding d which are 

d = i?(0)(M; 1 (s)p*) 12 - R(6)(A( s)s) 12 (3.58) 

Also for these two equations it is clear that there are two terms. One is always the geometric one, 
depending only on the shape s and s. The other one integrated over a gives 


/ -R(0O))( M r ( Q! ( T ))P*( T ))i 2 dr ( 3 - 59 ) 

Jo 

which is due to the presence of the impulse. 

The two additional terms (3.571 and (3.591 are exactly the so called dynamic phase presented in 
section 2.2 and represent the gap on the fiber (di,d 2 , 0 ) performed by the swimmer after a periodical 
change of shape. 
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4. Controllability 

In this section we will focus on the controllability of our system, i.e. its ability to move everywhere in the 
plane changing its shape. 

First of all we introduce some classical definition and results that will be useful in what follows 

4.1. Tools in geometric control theory 

Let us consider the following control system 


q = F{q,u) 


(4.1) 


where q are local coordinates for smooth manifold Q with dimQ = n and u : [0,T] —► U C M m is the set of 
admissible controls. The unique solution of (4.1) at time t > t 0 with initial condition q(t 0 ) = q 0 and input 
function it(-) is denoted q(t, to, qo, u). 

Definition 4.1. • The reachable set 1?' (qo,T) is the set of points in Q which are reachable from qo at 

exactly time T > 0, following trajectories which, for t <T remain in a neighborhood V of qo 

• The system is locally accessible from Xo if, for any neighborhood V of qo and all T > 0 the set 
R)£(qo) = Ut<r-^ 1 (<7ojt) contains a non empty open set. 

• The system \4-ty is locally strong accessible from qo if for any neighborhood V of qo and all T > 0 
sufficiently small, the set R v (qo,T) contains a non empty open set. 

• The system \4-ty is controllable, if for every q \, qi £ Q exists a finite time T > 0 and an admissible 
control u : [0, T] —> U such that q(T , 0, qi,u) = <72 


Let now suppose the system (4.1) to be an affine non linear control system, namely 


q = F{q,u) = f{q) + ^gj{q)i 


(4.2) 


1=1 


We now present some general results for this type of control systems 

Definition 4.2. The strong accessibility algebra Co is the smallest subalgebra of the Lie algebra of smooth 
vector fields on M containing the control vector fields g\ .. .g m , which is invariant under the drift vector field 
f, that is [f,X] G Co, VX £ Cq, every element of the algebra Co is a linear combination of repeated Lie brackets 
of the form [X k , [X k -i, [■.., [Xi ,gj].. .]]] for j = l...m and where X { G {/, 51 ? • • • 5 9m} • 

The strong accessibility distribution Co is the corresponding involutive distribution Co{q) = {X(q)\X G 

Co}- 


Proposition 4.1. Let q e be an equilibrium point of the system (j.2). The linearization of the system (4-2) at 
q e is controllable if 


rank 


,<9/ I ^ ^ f ^2 


I \t d h n ~ 1 

s| -%) 9 


We say that the Strong Accessibility Rank Condition at qo G Q is satisfied if 

dim Co(qo) = n 


(4.3) 


(4.4) 


Proposition 4.2. We say that the system (4-2) is locally strong accessible from qo if the strong accessibility 
rank condition is satisfied. 


Proposition 4.3. If the system (4-2) is driftless, namely 

m 

q = T.Uigiiq) 


(4.5) 
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its controllability is equivalent to its strong accessibility. 


Let us recall the definition of iterated Lie brackets 15 


Definition 4.3. Let f £ C°° and g £ C°° we define by induction on k £N ad^g £ C° 


ad° f g := g 

ad) +1 g := [/, ad'jg], Mk £ N. 


We are now ready to give a sufficient condition for small time local controllability 

Theorem 4.4. Assume that the controlled vector fields gi ■ ■ ■ g m generate a Lie algebra Cie{gi---g m } that 
satisfies Zie{g\ ■ ■ ■ g m } = T q Q for all q in Q then the corresponding affine system 


m 

<? = /(<?) + ^9i{q)ui 


is strongly controllable whenever there are no restrictions on the size of the controls. 


4.2. Swimmer controllability 


Let us consider the control system (3.441, since they involve the impulse p* we have two different type of 


control system depending on the initial value of this impulse. If it is zero, we have a non linear drifltess affine 
control system, whose controllability can be proved with classical techniques, instead if it is not zero we have a 
non linear affine system with drift, which is more tricky to study. 

Remark 4.1 (Scallop Theorem). Note that in the case of zero initial impulse, if we have only one real shape 
parameter we are exactly in the case of the famous Scallop Theorem according to which if the swimmer 
performs a cyclical shape change a the net motion of the swimmer after a period is null. 


Ad N 

Ad. 


\ r T r a ( T ) 

I = / A(s(t))s(t) dt= A(a) da = 0 since a(0) = a(T) 

J JO Ja(0) 


(4.6) 


Now let us study the controllability of this system in both cases of interest: Pg = 0 and Pg 7 ^ 0. 


4.2.1. Case 


P*(0) = 0 


In this subsection we want to study the controllability of the system which starts with zero impulse. Accord¬ 
ing to what said before this means that we deal with a non linear driftless affine control system. 


Case of 3 real shape parameters 


In this section we study exactly the case of three real controls, then we will generalized the results obtained to 
a larger number of parameters. More precisely, suppose that the deformation of our swimmer is governed by 
Si, S 2 , S 3 and according to [32] its shape is described in polar coordinates in the body frame by 

F(o , s) = 1 + e(si cos(2cr) + S 2 cos(3<r) + S 3 sin(3u)) (4.7) 

The perfect irrotational fluid has density p and the potential if* can be determined solving the Laplace problem 
with Neumann boundary conditions following the steps described the preceding sections. 

After that it is possible to compute the expression of the connection and the equation of motion 
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/ d :\ 

e 


pt 


p 2 


P 3 


Si 


S 2 

\s 3 / 



-(1 -/i)s 3 
0 
0 
0 

-(1 -aO(s3Pi +S 2 P 2 ) 

1 

0 

V 0 


2 

e iti + 


+ 


/ ~ Sl \ 


/ ° \ 

0 


-Si 

27rps 3 


27rps2 

M 


M 

27rps3p 2 


2-KpS2V2 

M 


M 

2irps 3 p* 

£ #2 + 

2irps 2 p* 

M 


M 

S 1 P 2 


-S 1 P 2 

0 


0 

1 


0 


(4.8) 


2 

€ 


with p = and M the mass of our body. 


Due to the change of variables (3.8), the equations of motion have to be supplemented with a so-called 


reconstruction equation allowing to recover d knowing 9: 


/ d \ 


61 


Pi 


P 2 


P3 


Si 


S 2 


\s 3 / 



\ 


9 Af 

27rps3p 2 

M ' 

27rps3p^ 

M 

S1P2 

0 

1 

0 


2 

e U 2 + 




(4.9) 


-si 

_ 27TpS2 

x Af ' 
27rps2P2 

M * 
27rps2P^ 

M 

-SlP 2 

0 

0 

1 


2 

e M 3 


Theorem 4.5. TTie system (4-9) is controllable. 


Proof: First of all note that system (4.91 is clearly of the type 

3 


q = '^2g i (0,p*,s)u i 


i =1 


Since the initial impulses are zero it is reduced to only six non trivial equations, indeed we easily have that 

p *(t) = 0 Mt 


































30 


TITLE WILL BE SET BY THE PUBLISHER 


(3.541. 


is a solution of the equations regarding p* 

Accordingly to theorem (fOl to prove the controllability it suffices to verify the Lie algebra rank condition, i.e 


dim(Cie{gi}i = 1,2,3) = 6. We compute all the vector fields gi and the Lie brackets of the first order [gt,gj\ with 
i ^ j (details in the Appendix) and compute their determinant 


det|gi,ff 2 , 33 , [91,92], [ 92 ,93], [51,53] | 


47 igpe 18 (gM - 2n(g - 1 )p (si + s|)) 


M 2 


(4.10) 


which is not null except for values S 2 = S 3 = 0. Since we can control the shape parameters we are always able 
to move from these configurations, and therefore cross the submaifolds defined by the equations S 2 = S 3 = 0 . 
Thus we can conclude that 51 , 52 , 53 , [ 51 , 52 ], [ 52 , 53 ], [ 51 , 53 ] are always linearly independent and dim(Lie{gi,i = 
1,2,3}) = 6, which proves the controllability result. □ 

General Case: m > 3 

In this subsection we deal with a generalization of the previous controllability result. Suppose that the shape 
of the swimmer is described by m real parameters Sj, i = 1 • • • m, which define a transformation near to the 


identity, whose expression is a generalization of formula (4.71. Moreover recall that we are still in the assumption 


that the swimmer starts with zero initial impulse in body coordinates. In this case the equation of motion turn 
out to be 

/ x \ 

y 
9 


Si 


\Sm ) 


= 2^ 9^ 

i =1 


(4.11) 


Note that also in this case, since the initial value of p* is null, p*(f) = 0 is still a solution and therefore p* does 
not appear in the system. We now investigate the controllability of the system (4.111. 


Theorem 4.6. The system (4-H'l is controllable. 


Proof: First of all observe that if we keep constant and equal to zero the last m — 3 controls, i.e. u t = 0, 7 = 
4-- -to the last m — 3 equations gives us easily s, (i) = 0 Vf, Vto > 4. This means that the shape of the 
swimmer is actually described by only 3 parameters. Therefore the remaining control equations have to be the 
same of the ones obtained in the previous section with m = 3. This implies that the first six components of the 
vectors gj\ Si =oi= 4 --m, j = 1,2,3 have to be equal to the vectors gi defined before. As a consequence we have 
that 


Lie{ ( g 1 ) , i = 1, 2, 3} C Lie{gt, i = 1 • • • m} 


(4.12) 


Moreover we have also that the vector space generated by the last m — 3 vector fields gi evaluated at s, = 0, i = 
1,2,3 have to be contained in the Lie algebra generated by all the 5 ^, since they are some of the generators. 


span{gj\ Si = 0 ,i=i,2,3, j > 4} C Lie{gi, i = l---m} 


(4.13) 


Furthermore we have also obviously that 


Lie{ 



1,2,3} n span{gj\ Si = 0 ,i=i, 2 , 3 , j > 4} = {0} 


(4.14) 
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This implies 


dim(^Lie{gi , i = 1 • • • m}^ > dim(^Lie{ , * = 1,2,3}^j + 

(4.15) 


>ra—3 

where the first equality derives from the proof done before in the case m = 3. 
Thus finally we obtain that 


=6 


dim (span{cjj | Si =o, i=i, 2 , 3, j > 4}) 


dim 




>m + 3 


(4.16) 


which proves the controllability of the system. 


□ 


4.2.2. Case pg ^ 0 

Let us suppose that our deformable body has an initial constant impulse Pg that is not null. As a consequence 
our control system is a system with drift of dimension m + 6. 

Case of 3 shape parameters 

We start we the simplest case of three control shape parameters; Since we start with an initial impulse that is 
not null we have the following control system with drift 


A 


/ M r (s)p* \ 

6 
• * 


(M7 1 (s)p*) 3 P2 

pi 
• % 
p 2 


(M r 1 (s)p*) 3 p* 

P 3 

— 

(M“ 1 (s)p*) 2 pt - (M~ 1 (s)p*)ip 2 

s\ 


0 

S2 


0 


\S 3 


J 


+ 


( ~{l~g)s2 

-(1 - g)s 3 

0 

0 

0 

-(1 - g)(s 3 pl + s 2 p* 2 ) 

1 

0 

V 0 


\ 


/ 


Z i 

e u i + 


-si 

0 

27rps 3 
M 

27rps3p 2 

M ' 

27rps3p 1 

M 

SlP2 

0 

1 

0 


\ 


2 

e U 2 + 


/ 0 \ 

-Sl 
_ 27TPS2 
M 

2-KpS2P2 

M * 
27rps2Pi 

M 

-S1P2 

o 
o 
i 


V 


(4.17) 


Which taking into account the reconstruction equations becomes 
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/ d \ 


e 


i>i 


P 2 


t >3 


Si 


S2 

. 



7l(6>)M7 1 (s)p* 

(M r 1 (s)p*) 3 p 2 

(M r 1 (s)p*) 3 pi 
--1 


(s)p *)2 P *1 - 
0 
0 
0 


(s)p*)iP2 


+ 



(4.18) 


Theorem 4.7. The system ( 4 . 18) is strongly controllable if there are no restrictions on the size of the controls, 
except at least on submanifolds of co-dimension greater than one defined by p\ = 0, p% = 0, p 3 = constant ^ 0. 
Moreover these submanifolds are invariant and the control system restricted to them is strongly controllable if 
there are no restrictions on the size of the controls. 


Proof: The system (4.18) is clearly of the type 


3 

9 = f('?) +^2&i(d)ui 

i=1 


(4.19) 


Applying theorem (4.4) to prove the strong controllability, we have to verify that the Lie algebra generated by 
the vector fields g 4 has the same dimension of the tangent space, i.e dim(Lie{gi,i = 1,2,3}) = 9. Thus we 
compute the Lie brackets of zero, first and second order of the vectors g^ (the detailed expressions are in the 
Appendix). 

The determinant of these vector fields is 





















TITLE WILL BE SET BY THE PUBLISHER 


33 


det jgl, g 2 , g3, [gl, g 2 ], [gl, gs], [g2, gs], [gl [g2, gs]], [g2 [g2, ga]], [g 3 [g2 , gs]] | 


8192 7 , 7 2 2 36 /,,r * n * 

7T ftp s 2 s 3 e ( Mp 2 - 2np 1 ps 2 s 3 ) 

(^M 2 (pl((2(p - 3)p + 3)s 2 - ps 3 ) + +p 2 ((2p - 3)s 2 + ((9 - Ap)p - 6)s 3 ))- (4.20) 

- 2n(p - l)Mp(2ppls 2 (s| - 2sl) - pl(s 2 + s 3 ) (4s| - 3s 2 s 3 + s 3 ) - 

- P 2 S 3 (-2 psl + si + sj)) + 8n 2 (p - l) 2 p 2 s 2 s 3 (s| + s 3 ) {jp 2 s 2 - p*s 3 )) • 


This is not vanishing except the following cases 






s 2 = 0 

This case is easily solved, indeed we are controlling s 2 , therefore we are always able to move from this 
configuration and cross the hypersurface s 2 = 0 . 
s 3 = 0 

This case is solved exactly in the same way as the preceding one. 

P * 2 = 0 

We want to prove that we always have a vector field that is non tangent to this hyper surface. For 
this purpose we compute the scalar product between the gradient of the determinant and each vector 
field g,; and see if it is non zero along the hyper surface p 2 = 0. Let us consider for example the scalar 
product with g 2 


((Vdet{gi,---},g 2 )|p* = o)M n = 

32768?r 9 p(pt) 3 p 9 S 2 S 3 e 38 ^ M 2 ((—2p 2 + 6p — 3) s 2 + ps 3 ^ + 2nMp(2 (p 2 — 3p, + 2) s 2 + 
+£ 2(^3 — ps 3 ) — 2 (2p — 3p + l) s 2 s 3 — (ji — l)s 3 ) + 871 (p — 1 p s 2 s 3 + 53 )^ 


this 

a) 

b) 

c) 


is null only in the following cases 
s 2 = 0 


From which we are always able to move as we have seen before 
s 3 = 0 

From which we are always able to move as we have seen before 

V\= 0 

Looking at equations (4.171 this case means that we are on the submanifold p\ = 0, p 2 = 0, 
p 3 = const that has co-dimension 3. In this case the system cannot leave this submanifold, which 
is invariant. 


d) s 2 = f (s 3 ) 

Here / is a suitable function of S 3 such that 

{^M 2 ( (—2 p 2 + 6p — 3) s 2 + ps 3 ) + 2nMp{2 ( p 2 — 3p + 2) s 2 + 

+ s 2 (s 3 - ps 3 ) - 2 (2 p 2 -3/i + l) s 2 s 3 - {p - l)s 3 ) + 8n 2 (p - l) 2 p 2 s 2 s 3 (s 2 + s 3 )^) = 0 


Since we are controlling both s 2 and s 3 we can always move from this configuration. 

In any case the points in which this scalar product is zero define a submanifold of co-dimension grater 
than one that we are always able to bypass. Therefore we have proved that the vector field g 2 is non 
tangent to the hyper surface p 2 = 0, except in the case of p* =0, thus we are able to move from it using 
suitable controls. 

» _ 2t Tp\ps 2 s 3 


Like in the previous case we compute the scalar product of the gradient of the determinant with g 2 and 
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see that it is non null except on sub manifold of co-dimension grater than one, more precisely 
«Vdet,g 2 >| . 2 „ p> s 2 s 3 )M 14 = 

p 2 =- 1U - 

655367r 9 /.ipi p 9 s|s3e 38 (M 2 + 27 r 2 p 2 S 2 S 3 ) (A/ 3 (( 2 p 2 — 6p + 3 ) S2 — /.1S3) — 

2 nM 2 p (2 (p 2 — 3/i + 2) s 3 + (4 — 3 p)« 2 S 3 + (4 — 3 p)s 2 S 3 — (p — l)ss) — 

47r 2 (p - l)Mp 2 S2sl ((4p - 3) s 2 2 + (2p - 3)s 2 ) + 167r 3 (p - l) 2 p 3 s|s 2 (s| + S 3 )) 




this is null only in the following cases 

0 S 2 = 0, S 3 = 0, p\ = 0, S 2 = h(s 3 ) where h. is a suitable function of S 3 in which the scalar product 
vanishes. 

These cases have been already faced and treated before. 

P *2 = f(Pl,S 2 ,S 3 ) 

Where / is a function that represents the value of p 2 in which the last factor of the determinant is 
null. Also in this case computing (Vdet,g 2 ) | p » = j( p » S2 S3 j we see that all the points in which it is null 
define sub manifolds of co-dimension grater than one defined by the equation p\ = 0, p 2 = 0 and 
p 3 = const ^ 0 . 


Finally we can conclude that the vector fields 


[§1) § 2 ]; [Si; S 3 ] 7 [f?2 7 S 3 ] 7 [Si [S2 7 S 3 ]] 7 [S2 [S2 7 S 3 ]] 7 [S3 [S2 7 S 3 ]] 


are always linearly independent except the submanifolds defined by pi 2 = 0, p\ = const, which proves the first 
part of the theorem. 

These submanifolds are invariant and on them the dimension of the system reduces to 6. If we restrict to 
them, we have still a system with drift, whose strong controllability can be proved verifying that the vectors 
<717 <?27 <?37 [gi,g 2 ], [<727 33]7 [<7i 7 33 ] are linearly independent (see Theorem 4.4). Therefore the system restricted to 
each of these invariant submanifolds is strongly controllable if there is no restrictions on the size of the controls. 

□ 


General case m > 3 


In the case of initial impulse not zero, as we have said before we have a control affine system with drift of 
dimension m + 6. 

fdt\ 

d * 2 

0 

Pi 

P *2 =f + V)g im (4.21) 

P*3 i=1 

Si 

\Sm/ 


Theorem 4.8. The system (4-21) is strongly controllable except at least on submanifolds of co-dimension greater 


than one defined by p* = 0, p 2 = 0, p$ = constant 7 ^ 0. These submanifolds are invariant and the control system 
restricted to them is strongly controllable if there are no restrictions on the size of the controls. 


Proof: To prove the strong controllability of the system (4.211 we exploit the theorem (4.41. Thus the condition 
to prove is that the dimension of the Lie algebra genereted by the control vector fields g; has dimension m + 6. 
The proof is similar to the proof of the previous theorem. Using the following facts: 
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• Lie{ ,i= 1,2,3} C Lie{g t ,i = 1,- • ■ ,m} 

• > 4} C Lie{g h i = l,--- ,m} 

• Lie i (o) ,* = 1,2,3} n span{gj | « 4 =o, *= 1 , 2 , 3 , j > 4} = {0} 
we deduce that 

dim(^Lie{gi, i = 1 • ■ ■ to}^ > dim(jLie{ , i = 1,2, 3}^j + 

(4.22) 


>m—3 


=9 


+ dim^span{gj| g4 =o i i=i j 2,3, j > 4}) 


Which proves that dim^Lie{gj, i = 1 • • • m}J = m + 6, except at least on the same submanifold on which the 
condition is not satisfied even in the case m = 3 that are of co-dimension greater than one defined by p\ = 0, 
P 2 = 0, P 3 = constant yf 0. Again if we restrict ourselves on this invariant submanifolds, recalling that in 
Theorem |4.6| we have proved that the dimension of the Lie algebra generated by the g; is to + 3, the condition 
of controllability of Theorem |4.4| is satisfied. Therefore the system is strongly controllable on the invariant 
submanifolds if there are no restrictions on the size of the controls. 


□ 


Conclusions and perspectives 

In this paper we have investigated the geometric nature of the swimming problem of a 2-dimensional de¬ 
formable body immersed in an ideal irrotational fluid. 

We faced a new problem: the study of the controllability properties of a dynamical system which can start 
with a non zero initial impulse. Reinterpreting the hydrodynamic forces exerted by the fluid on the body, as 
kinetic terms, and describing the shape changes with a finite number of parameters, we derive the equation of 
motion of the system. Using classical techniques in control theory we are able to gain some good results for the 
controllability of this kind of system. 

If it starts with zero initial impulse we recover results present in the literature. We are always able to find a 
suitable rate of deformation which makes the swimmer moving between two different fixed configurations. If 
instead the body starts with an initial impulse different from zero, the swimmer can self-propel in almost any 
direction if it can undergo shape changes without any bound on their velocity. 

The fact that we take into account the presence of an initial impulse not null, and the analysis of the control¬ 
lability of this system seems innovative and makes the study of the self-propulsion of deformable bodies in an 
ideal fluid more accurate and complete. 

The approach described in this paper can be extended in a number of natural ways. To begin with, we have 
restricted our attention to planar swimmers. The general 3-dimensional case is conceptually straightforward, 
even though the way of describing the shape changes should be different. 

The study of bodies that change their shape using only a finite number of parameters is the initial point of 
a more complex study of controlling the deformation by diffeomorphisms. Future work will also explore the 
optimal control problem associated to these kind of systems, especially in the case of non zero initial impulse. 
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Appendix 


The vector Fields gi and their Lie brackets of the first order mentioned in theorem 4.5 


are 


3i 


9 3 


1Z(6) 


V 

/ 


no) 




/ 

/ -Si 

\\ 



o 




\ 27TPS3 

/ 

92 = t 


\ M 

0 




1 



V 

0 

/ 


(4.23) 


The Lie brackets generated by these vector fields are 


[31,32] 


[31,33] 


[ 32 , 33 ] 


/ 


V 

/ 


V 

/ 


2-7TS2S3 (/A — l)p sin(0) — cos(0) (M/A — 27TS3 (/A— l)p) \ 

M 

sin($) ( 2-7733 Cm - l)p—A-f/ a) — 27732^3 (/a— l)p cos(d) 

M 

0 

0 

0 

0 / 

sin(@) (Mp — 2-7732 (/A— l)p) — 2-7732 33 (/A — 1)P COS(0) \ 

M 

cos(0) ( M/A — 2-7732 (m — 1)c)+ 277S2 33(/A— l)p sin($) 

M 

0 

0 

0 

0 / 

27rsip(s2 sin(0)+S3 cos(0)) 

M 

27rsip(s3 sin(0) —S 2 cos(0)) 

¥ 

_ 47rp 

M 

0 

0 
0 


(4.24) 
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The vector fields that we need to compute the Lie algebra generated by g; in theorem |4.7|are 


gi = e 


g3 = e 




g2 = e 


■2TTpS 3 p% 

M 

2irps 3 pl 

M 

SlP2 

0 

1 

0 


(4.25) 
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Their Lie brackets of the first order are 


[gl,g2] 


[gl, gs] 


[g2,g3] 


/ 27TS2S3 (//, —l)p sin(0) —cos(0) 27TSg (/i —l)p^ \ 

' M ' 

sin(0) ^27TSg (/i — l)p — M — 27TS2S^(fx — l)p cos (6) 

M 

0 

0 
0 

27TS3 (// —l)p(p^S2 —P 2 S s) — Mp2 (m — 2) 

M 

0 

0 
0 

/ sin(0) —27TS2 (/u. —l)p) —27TS2S3 (/^ —l)p cos(0) \ 

I M ' 

cos(0) ( M p — 27TS2 (m — 1)p) +27TS2S3 (£t — l)p sin(0) 

M 

0 
0 
0 

(27rs 2 (e- 1 )p(P*S2+p5 s 3)--iH(p*(^-l)+P2)) 

M 

0 
0 
0 

/ 27TSip(s2 sin(0)+S3 cos(0)) \ 

/ M \ 

27rsip(s3 sin(0) —S 2 cos(0)) 

¥ 

47rp 

AT 

47rp^Mp2 —27r Pi s 2 s 3P) 

77 ^ ; 

87T p 2 S2 s 3P 

7 A* 2 

27rp 1 sip(s2+S3) 

M 

0 
0 
0 


(4.26) 
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Finally the only non zero brackets of the second order are 


[gl, [g2,g3]] 


[g2, [g21 g3 ] ] 


[g3 ■ [g21 g3]] 


( 


6 


6 


V 

/ 


e 


6 


V 

( 


e 


6 


V 


27r(2/x —3)p(s2 sin(0)+S3 cos(0)) 

M 

2ir(2fi — 3 )p(s 2 cos( 0 ) —S 3 sin( 0 )) 

M 

0 

0 

0 

27rp(M(pJ(s2+S3) + 2p2 s 3(M—l))+4irS2S3(M—1)P(P2 S 2—Pl s 3)) 

IT 2 

0 

0 

0 


27rsip^sin((?)^3M—27TSgp^ + 27TS2P3P cos( 0 )) 



M' 2 

0 


I 67 r 2 s 3 p 2 ( M p \ + 277^2 S2«3p) 

T- 

167T 2 s 3 p 2 ( Mp* —2irp*S2S 3 p) 

i ^ ; 

27rs 1 p^Mp 1 +27rp 2 s 3 p(3s2+S3)J 
U 2 

0 

0 

0 

27rsip(cos(P) ^3 AT—27rs|p^ +27TS2 S3 p sin(P)^ 
M 2 

2ttsip( 3M sin(0) —27TS2 p(s 2 sin(0)+S3 cos(0))) 

P 

0 


167T 2 S2p 2 [MpJ —27rp2®2S3P^ 

IP 

167r 2 S2p 2 (Afp 2 — 27rp* S2S3p^ 
IP 

27rsip(Mp* — 27TP2S2P( S 2+3S3)^ 

IP 

0 

0 

0 


\ 


/ 

\ 


) 


\ 


) 


(4.27) 
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